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Abstract. The Fisher-Bingham system is a system of linear partial differ- 
ential equations satisfied by the Fisher-Bingham integral for the n-dimensional 
sphere S n . The system is given in [4j Theorem 2] and it is shown that it is a 
holonomic system [T] . We show that the holonomic rank of the system is equal 
to 2n + 2. 
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1 Introduction 

Let x = (xij) and y = (yi) be parameters such that xtj = Xji for i ^ j. Let 
Z be a function, which is the normalization constant of the Fisher-Bingham 
distribution, defined as 

n+l 

Z(x,y,r)=l exp | y] XjjUtj + y^ y Viti \ \dt\ (1) 

i=i 

where S n (r) = {(t\, . . . , t n +\) \ Y^i=x *i = r2 ' r > 0} is the n-dimensional 
sphere and \dt\ denotes the Haar measure on the sphere. 
Let D be the Weyl algebra 

D = C( 

where dij = d/dxij, dk — d/dyk and d r = 3/ dr. It is shown in |T] and j4] 
that the normalization constant (fT]) of the Fisher-Bingham distribution is a 
holonomic function in x, y, r and consequently it is annihilated by the following 
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holonomic ideal I in D generated by the following operators in D: 

dij - didj, 

n+l 

i=i 

+ 22 (xsjdid s - Xisdjds) + yjdi - yidj, 
rd r - 2 E x^didj - 2J vA - n. 

i<j i 

We call the system of differential equations defined by / the Fisher- Bingham 
system. 

For a left ideal J in D, the holonomic rank of J is defined as the dimension 
of the K = C(xij ,yk,r | 1 < i < j < n + 1,1 < k < n + 1) vector space 
K{dij,dk,d r )/K(dij,dk,d r )J. The rank is denoted by rank (J). When J is 
a holonomic ideal, the rank is finite. The holonomic rank agrees with the di- 
mension of the holomorphic solutions of the associated system of linear partial 
differential equations at generic points and with the size of the Pfaffian equation 
associated to J. As to general facts on the holonomic rank, we refer to, e.g., 
the chapters 1 and 2 of [6 . The holonomic rank is a fundamental invariant 
of the D- module D/J and there are several attractive studies on holonomic 
ranks. For example, Miller, Matusevich and Walther studied holonomic ranks 
of A-hypergeometric systems by introducing a new homological method [3] . 

We are interested in the holonomic rank of the Fisher- Bingham system /. 
We prove the following theorem in this paper. 

Theorem 1 

rank (I) = In + 2 

In [3] , we proposed a new method in the statistical inference which is called 
the holonomic gradient descent. The method utilizes a holonomic system of 
linear partial differential equations associated to the normalization constant. 
The complexity of the method depends on the holonomic rank and correctness 
of the method are proved by utilizing the holonomic rank. In the case of the 
Fisher-Bingham distribution, which is the most fundamental distribution in the 
directional statistics, Theorem [T] is applied in [2], which gives a generalization 
of the result in [3] shown with a help of a computer program. Our method 
to prove the theorem is the Grobner deformation to the direction (— w, w), 
which is discussed in [B] for A-hypergeometric systems, and a determination of 
Grobner bases by hand with adding several slack variables which do not change 
the holonomic rank. 
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2 The Rank of the Diagonal System 

When the matrix x is diagonal, the normalization constant Z satisfies a system 
of linear partial differential equations for the variables xu, yk, r. Let / be the 
left ideal in D generated by 

A l = d u -Of (l<i<n + 1), 

n+l 

Cij = 2{xu - .i- ,, )(),(), + y t dj - yjdi (1 < i < j < n + 1), 

n+l n+l 

E = rd r - 2 ^2 x a d i - Vidi ~ n 

i=l i=l 

and dij, i ^ j. The ideal I annihilates the function Z restricted to the diagonal 
of x 0. 

Theorem 2 The holonomic rank of I is 2n + 2 . 

Our proof of the theorem reduces to the proof of the following proposition. 

Proposition 1 Let R be the ring of differential operators with rational function 
coefficients 

R = C(xn, x n+ln+ i,y!, y n+1 ,r)(d n ,. . . , d n+ln+ x, d 1} d n+ i,d r ). 

Let RI be the left ideal of R generated by A4, B, Cij, E. Let < be the term order 
on R which is the block order with d r 3> {da} 3> {dj}. The order of the block 
{da} is the graded lexicographic order with du > • • > <9„+i n +i and that of the 
block {di} is the graded lexicographic order with d\ > • • • > d n +\. A Grobner 
basis of RI with respect to the term order < is 

Ai = da - df (i = l,...,n + l), 

n+l 

i=i 

= 2(x it - Xjj)didj + yidj - yjd t (1 < i < j < n + 1), 
( We put a,, = 2(x u - Xjj),Fij = y.i), - yjdi,Cij = n,,0,i), + F^), 
D k = d k B - dia'^Cik d k -!a^ lk C k - lk (k = 1, . . . , n + 1), 

n+l n+l 

E = rd r - 2 ^ xudf - yidi - n. 

i=l i=l 

The initial monomials of the Grobner basis are 

in<(Ai) = in< in< (5) = in < (ai) 2 ,in < (C y ) = in < (a t )in<(9 i ), 
in<(D fc ) = in < (a fc ) 3 ,in < (£;) = in<(9 r ). 
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Here, the initial monomial in<(c(x, y, r)dj L] d^" L] dk) is denned as the el- 
ement c(x, y, r)t 7 Y[ II Vk k m the polynomial ring with rational function 
coefficients C(a;, y, r)[fn, . . . , £ n+ i n +i, r/i, . . . , rj n+1) t] (see, e.g., Chapter 1]). 

By the proposition, the standard monomials of the quotient ring R/RI are 
1, dx, 82, 8f, • • • 1 9 n +i, dn+i- Therefore, the holonomic rank of I is 2n + 2 (The- 
orem [5]). Let us prove the proposition. 

Since D k is expressed by B, Cik, ■ ■ ■ , Cfc-ife, the operator D k is the ele- 
ment in RI. In order to prove Proposition [TJ we will show that any 5-pair for 
Ai,B, dj,Dk, E is reduced to by Ai, B, dj,Dk,E. The following lemmas are 
proved by straight forward calculations. 

Lemma 1 Let P and Q be elements in R. If the initial monomials are coprime, 

i.e., gcd(in<(P), in<((5)) = 1, then the S-pair S{P,Q) is reduced to [P,Q] by 

P and Q (we denote the reduction by S(P,Q) >* [P, Q]), where [P,Q] is the 

P,Q 

commutator of P and Q. In particular, when [P,Q] — 0, the S-pair S{P,Q) is 
reduced to 0. 

Lemma 2 We have 

[A p , Cij] = 0, 
[B, <?«] = (), 

[Cij , Cjk] = Cik, [Cij, Cik] = ~Cjk, [Cik,Cjk] — ~Cij [i < j < k), 
[C l3 ,C pq }=0 ({i,j}n{p,q} = <b). 

Lemma 3 We have 

r (» < j) 

[Di,Aj] = } 2aJ*djCji (i > j), 

lEti2o^C„ (i = j) 

[Di,B]=0, 

[A, A] = -Bdj + J2 a^diidiQj + ftCy) + a^o^^CJy). 

Z<i Ki 

When i,j ^ fc, we obtain 

(0 (k-Ki) 
[Cij,D k ] = lo (j<k-l) . 

{ - a 7k Py , d l Cik] = a^ 1 {d t C ]k +djC ik ) (i < k - 1< j) 

Lemma 4 We have 

[A h E] =0, 
[B, E] = -2B, 
[C lJ ,E]=Q, 

i-l 

[A, E] = -3A - 2 ^ H^d k C ki . 

fe=i 
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Proof of Proposition^ We prove that any S-pair for A i: B, Cij,Dk,E is reduced 
toO by Ai,B,C ih D k ,E. 

S -pairs of Ai and Aj, B,Cij. The initial monomials are coprime, and the 
elements commute. By Lemma [TJ we obtain 

S(Ai, Aj) — ►* 0, 
S{A U B) — ►* 0, 
S(Ai, Cjk) — ►* 0. 

S-pair of B and Cij. When i > 1, the initial monomials in<(i?) = in<(<9i) 2 , 
in<(Cy) = in<(9i)in < (9j) are coprime. Operators B and Cij commute by 
Lemma [2] By Lemma [TJ we have S(B,Cij) — >* 0. 

Wheni = 1, the initial monomials in< (B) = in<(9i) 2 , m < (Cij) = in < (9i)in < (9 3 ) 
are not coprime. We obtain the following reduction sequence of the S-pair: 

S(B, Cij) = aijdjB - diCij 

= a:MU)i + ■ ■ ■ + $ + • • • + dX+i fyr 2 ) d-\ i\ , 

= axjijdjdl + ■ ■ ■ + % + ■ ■ ■ + 8,81+1 - d,r 2 ) - dia^Fu) 



— ->* a^ddjdi + ... + % + ... + djdi +l - d^) - d ia ^F l3 - d 2 a~ F 23 ) 
023 

>* >* a Xj Dj — > 0. 

S-pair of Cij and Cm ({i,j}Ci{k, 1} — 0). The initial monomials in<(CV/) = 
in<(9 i )in < (9j), in<(Cfc;) = in < (9fc)in < (9;) are coprime. Operators CV,- and Cu 
commute by Lemma [21 By Lemma [TJ we obtain 

SiC^Cu) 0. 

S-pair of Cij and Cjk (i < j < k). We have the following reduction sequence 
of the S-pair: 

S(Cij,Cjk) — ajkdkCij — aijdiCjk = —aikyjdidk + ajkyidjdk + a^ytdidj 
- ~ ~ > * Vii-Fjk) + Vki-Fij) - yj{-F ik ) = 0. 

The S-pair of Cij and Cik and that of Cik and Cjk are also reduced to 0. 

S-pair of Di and Aj. The initial monomials in<(_Di) = in<(<9i) 3 , in.<(Aj) = 
in<((?jj) are coprime. When i < j, operators Di and Aj commute. By Lemma 
[TJ we have 

S(D u Aj) >* 0. 

When i > j, by Lemmas [TJ and [3j we have 

S(D i: Aj) >* [Di, Aj] = 2,, ~i),C,, >* 0. 

D t ,Aj Cji 
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When i = j, by Lemmas [T] and [3J we have 

i-1 

S(Di,Ai) — — ^* [A,^i] = ^2a i T 2 9 i a i — ^* 0. 

(=1 

S-pair of D, and B. When i = 1, the S-pair is S(D 1 ,B) = d x B - d x B = 0. 
When i > 1, the initial monomials in < (D i ) = in < (9i) 3 , in<(i?) = in < (9i) 2 
are coprime. By Lemmas [1] and [3J we have 

S(Di,B) >* [D i ,B}=0. 

S-pair of Di and Dj. The initial monomials in<(Dj) = in<(<9i) 3 , in<(I?j) = 
in<(9j) 3 are coprime. By Lemmas Q] and [3l we have 

s(D it Dj)^ n.n: ^>*o. 

S-pair of Cij and D k . The initial monomials are in<(Cij) = in < (9i)in < (9 :; ), 
in<(Z)fc) = in<(<9fc) 3 . When i ^ fe and j ^ fc, the initial monomials are coprime. 
By Lemmas [1] and [3l we have 

S(Cij,D k ) — — +* [C ij} D k ] — — ►* 0. 

When i = k, the initial monomials in<(Cy) = in<(9j)in<(S,), m < (Di) = 
in<(9i) 3 are not coprime. This case needs a care of an order of applying reduc- 
tions. We will reduce the (S-pair by Dj and then reduce remainders by C^'s. 

S{C ij ,D i ) = d 2 i C ij -a ij d j D i 

71+1 i-1 

= dfF tJ - n,,i),ir - aijdj ( did f - d ^ 2 - H ^Fii) 

i=i+l,/^J 1=1 

n+1 j-1 

dfFij + aiM J2 d i 9 i - d ^ ~ E fajfti) 

J 1=3+1 1=1 

n+1 i-1 

i=i+i,i^j i=i 
= -aijdi(di + ia^ +lj Ci + ij H h dj-iaJ^Cj-ij) 

i-1 i-1 

- a-ijdiiYl a ij lFl ^ + a v d i(Yl d i a u lF i*) 
i=i i=i 

i-1 i-1 

>* 



-Oijdi(%2 a h x Fij) + aijd^Y^ dia^Fn) 



1 = 1 1 = 1 

i-1 
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Since a { - l didiFij — a^didjFu > * 0, the 5-pair S(Cij,Di) is reduced 

Cij ,Cn-,Cij 

to 0. 

When j = k, the initial monomials in<(Cij) = in<(9j)in<(9j),in<(Z)j) = 
in<(<9j) 3 are not coprime. This case also needs a care of an order of applying 
reductions. 

S{C ij ,D j ) = d 2 j C ij -a ij d i D j 

n+l j-1 

= Fijidf + 9?) - a % M djdf - d,r 2 - £ fto^Fy) 

i=j+l l=l,tyi 
n+l n+l j-1 

-+* F^- Y, df + r^-aMY didi-fyr 2 - E 

!=l,i/i,j J=j+i 1=1, tyi 

n+l n+l j — 1 

= f«(- £ 9 ") + r2 ( a iA d 3 + Fi) - oya«( Y d ^ E W 

l=l,l^i,j l=j+l 1=1, tyi 

n+l n+l j-1 

^* F a(- E 9 ?)-*iA( Y E fajfti) 

' J l=l.l^i,j 1=3+1 l=l,tyi 

n+l j-1 J-1 

= (Y df)(-"iAdj-Fij)-Fij Y d f + a iA Y fay-fti 

1=3+1 1=1, tyi 1=1, tyi 

j-1 j-1 

^*-Fij Y + Y ha^Fu 

' J l=l,tyi l=l,ljii 

3-1 

= Y dmji-a^diFij+a^diF^). 

l=l,tyi 



Since a^OiF,, + a^Fy >* 0, the S -pair S(Cij,Dj) is reduced to 

0. 

S-pair of E and A*. The initial monomials in<(F) = in<(<9 r ), in<(.Aj) = 
in<(9ji) are coprime. By Lemmas [1] and SI we have 

S(Ai,E) y* [Ai,E] = 0. 

Ai,E 

S-pair of E and B. The initial monomials in<(i?) = in<(9 r ), in<(_B) = 
in<(i9i) 2 are coprime. By Lemmas [I] and |4l we have 

S(B, E) >* \B,E] = -2B — ►* 0. 

B.E B 

S-pair of E and Cij. The initial monomials in<(i?) = in<(9 r ), in<(Cy) = 
in<(9j)in<(9j) are coprime. By Lemmas [T] and |4j we have 

S{Cij,E) — — >* [Cy,F]=0. 



S-pair of Di and E. The initial monomials in<(i?) = in< (d r ), in< (Di) = 
in<(9i) 3 are coprime. By Lemmas [T] and 21 we have 

i-l 

S(Di,E) — - >* [A,£] = -3A -2^0^^ — ; — -+* 0. 

fc=l 

We have proved that any S'-pair is reduced to 0. By Buchberger's criterion, 
the set {Ai, B, Cij, Dk, E} is a Grobner basis of RI. Q.E.D. 



3 Grobner Deformation of the Fisher-Bingham 
System 

Consider the system of differential equations I ■ f = 0. Intuitively speaking, 
we want to prove that the system / can be deformed to the diagonal system 
I without increasing the holonomic rank. This can be done by a Grobner ba- 
sis computation with a weight vector (—w,w) [6j Theorem 2.2.1]. However a 
straightforward calculation does not seem to be easy. We need to use some 
technical tricks to determine a suitable Grobner deformation. Since these tricks 
may look too technical for the general n, we explain them in the case of n = 1 
in Section @] to clarify our idea without technical details of this section. Readers 
are expected to refer to the Section 2] when technicalities get complicated. 

We will introduce new indeterminates to make our Grobner basis compu- 
tation possible by hand with employing the idea of the proof of [5J Theorem 
3.1.3]. Let a pq , hi, Cj, d (1 < p < q < n + 1, 1 < i < n + 1) be constants, 
which we call slack variables when they are regarded as indeterminates. We 

put g = yr d3 Yl p <q x pq q ^J f an d make a change of the independent variables 

Hi by yi + biCi. Then, the system of differential equations for the function g 
is I' • g = where I' is the left ideal in D generated by the set of operators 
G' = {A' pq ,B,C' l0 ,E'} where 

^PQ — X pq^pq ~ Xpqdpdq Qpq , 
n+1 

2=1 

+ ^ i x sjdid s - x is djd s ) + (yj + bjCj)d{ - [y l + kci)dj, 

n+1 

E' = rd r - 2 2J Xijdidj - }^(yi + hc^di -n- d 3 . 

i<j i=l 

The key fact is that the holonomic rank of / for / agrees with the holonomic 
rank of I' for g for any constants a pq , hi, Cj, d. 
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Let us make the same change of the variables for the diagonal system; let I' 
be the left ideal of D generated by the set of operators G' — {A' U ,B, C^E', Xijdij — 
a ij (* j)} where 

A' u = x H df - xudu -a% (1 < i < n + 1), 

n+l 

i=l 

Clj = 2(xjj - xu)didj + (yj + bjCj)di - (y t + hc^dj (1 < i < j < n + 1), 

n+l n+l 

E' = rd r -2^ xudf - J2(Vi + b i*) d i ~n-d\ 

i=l i=l 

The holonomic ranks of I and I' agree. 

Define the weight vector w by Wij = 1, (i ^ j), wu — 0, Wk — 0, and w r = 0. 
Here, stands for dij, Wk stands for dk, and w r stands for d r . The initial 
form inr W!W )(i), I S D, is the sum of the highest (— w, u>)-degree terms in I and 
m (-w,w){I') is the left ideal generated by I £ I' where the weight — w stands 
for space variables xy, yk, r corresponding to differential operators dij, dk, d r 
respectively (see, e.g., [BJ Chapter 1]). 

Theorem 3 For generic complex numbers a pq , hi, Ci, d, we have 

in(-«,,«;)(-0 = I'- ( 2 ) 

In order to prove the theorem, we regard a pq , bi, Ci, d as ring variables with the 
weight and consider the following homogenized system of 

— hXpqdpq Xpqdpdq ®pq, 
n+l 

i=l 

C^j - Xijd^ ~t~ ^i^Xjj Xn^jOiOj XijOj 

+ ^2 (x s jdid a - x is djd s ) + (hy 3 + bjCj)di - (hy l + biCi)dj, 

n+l 

E' h = hrd r -2^2 x ijdtdj - J2( h Vt + hci)di - nh 3 - d 3 . 

i<j i—1 

Let I' h be a left D h [a,b,c,d)-ide&\ generated by the set of operators G' h := 
{A'p q , B,C'tj , E' h }, where D h [a, b, c, d] is the homogenized Weyl algebra of the 
ring D[a, b, c, d] — C[a, b, c, d](xij, yk, r, dij, dk, d r ) with the homogenization vari- 
able h (see, e.g., Section 9]). We introduce a new term order <^_ w w \ over 
D h [a,b,c,d], which compares the total degree first, (— w, w, 0)-degree second, 
otherwise we apply the following block order as a tie breaker: d 3> r ^> {a pq \ 
i < j} > {b k } > {c fc } > {y k } > d r > {d^ \ i < j} > {d u } > {d k } > {x t] 
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i < j} 3> { x u} ^ h- Here, the block {bk} has a lexicographic order so that 
bi > 62 > • • • > b n+ \. The use of this tie breaker is a key of our calculation. 
Although, the initial monomial in^h (£) is an element of the associated 

commutative ring, we denote it by the associated element in D h [a, b, c, d] as 
long as no confusion arises. For example, we denote in < h (dij) by dij in- 
stead Of £ij . 

Proposition 2 A Grobner basis of I' h with respect to <^_ w w \ is G = 

We need four lemmas for proving the proposition. These can be obtained 
by a straightforward calculation. 

Lemma 5 The initial monomials of the generators of I' h are 

m< (- w u, o) (^pq) ~ ~ a pq' 

2. in^h (B) = -r 2 , 

3. ■m < H_ M (C$) = -b i c i d j , 
J. in.h (E' h ) = -d 3 . 

C — W,w0) 

In particular, they are pairwise coprime expect in the case that the pair in < h ^ (C^) 
and in < h 0) {^ik) and the pair (C'ik) and in ^ ^'Jk)- 

Lemma 6 The commutators of two generators of I' h are 

1. [A%,A%] = [A%B] = [A%,C%] = [A%E' h ] = (for any p, q, r, s, i, j), 

2. [B,C%] = 0, [B, E' h ] = -2hB, 

3- [C#,Cjg]=0 ({i,j}n{fc,/} = 0), 

; \nih f~flh] _ (~>ih . /~iih \r11h f~f>h] _ ruh \rnh (~ilh\ _ (~ith (a A ^ 1A 

4. [Oy , O jfe J — o fcl .— — O ife7 [Oy , L/ ifc J — l^ ]k , L^ifc ' ^rffeJ — U i] I* ^ J ^ K )> 



5. [C'^E' h ] =0. 



Lemma 7 T/ie following holds for S -pairs of generators of I th : 

1. S{A%, A%),S{A% B),S{A%, C$), S{A% E' h ) — ►* {for any p, q, r, s, 

2. S{B, C'lf) — >* ; S(B, E' h ) — >* 0, 



3. s(C#,C#)— >*o ({i,j}n{M} = 0), 

* S(C#,C#)— ►*() (t<j<fc), 



5. S(C%,E' h ) ^* 0. 



Proof. Lemma Q] holds in the homogenized Weyl algebra D h [a,b,c,d] too. 
Therefore, these follow from Lemmas [5] and [H Q.E.D. 
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Lemma 8 We put C$ := b jCj d t ~ b^dj and C% := C$ - Cfi. Then the 
following cyclic relations hold: 



1. 


dkC* + 




d 3 &t = c>!;d k + 




+ c' k % 


= 0, 


2. 


dkC% + 


be* + 


tfi'ik = C%d k + 




+ c' k h A 


= o ; 


3. 


dkC% + 


QiC% + 


djC'& = C%d k + 




+ C&% 


= o, 


1 




+ b i c i C% + b j c j C> k \ = C% 


bkCk + 







Proof of Proposition fj| By Lemma we only need to check that the S'-pairs 
S(C$, C'£) and S{C$, Cf k ) are reduced to zero. 

The former is S{C'^C'^) = d k C$ - djC% = d k C$ + djCfa, because the 
initial monomials arc in^h (C'/, L ) = —biCtd^ and m^h (C'A) = —bsCidtf 

The S-pair is equal to —diC'Jl by the [3j:d formula of Lemma [8l This implies 
that it is reduced to zero. 

The latter is S(C^, C'Jl) = bjCjCfy — &iC,C^, because the initial monomials 
are in^ ^ (C&) = -b iCl d k and in^ ^ = -b jCj d k . 

Firstly, we show that it can be expressed as 



s{c&, c%) = 53(«y fc> + l)x ks d s + hy k + b k c k c§ 



'n+l 



J2($is + l>isd s + h yi C% 



\s=l / 

+ (j2(S js + l)x js d s + hy^ Ct (3) 

Here, the symbol 5 is the Kronecker delta. This expression of the S'-pair is a 
key of our proof. Since each monomial of the left hand side (LHS in short) has 
at least one variable in 6^, bj or b k , we divide the right hand side (RHS in short) 
into two parts: Si contains these variables, S2 does not contain them. Then, 

/n+l \ /n+l \ 



+ ( + 1 + hyA C' k \ 
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is equal to zero by a straightforward calculation and we have 



St = b k c k C% + \J2( S ks + lWsd s + hy k j C, 

]{S ls + l)x ls d s + hyA C% + (j2(S js + l)x js d s + hyj \ C] 




\s=l 



l)x ks d s d l + hy k dij - ha \^_J$ ks + l)x ks d s dj + hy k d 3 j 

l)x is d s dj + hyidj^j - bjCj ^^(<5 ls + l)x is d s d k + hy^kj 

+ hci \^2(Sjs + l)x js d s d k + hyjd^j - b k c k (^2(S JS + l)x js d s di + hyjd^j 
— bjCjC'*}. - btCiC'j'l = S(C' ik , C'j'l). 

Finally, we show that the expression ([3]) is a standard representation. We 
note that the RHS of has three terms, which appear in the first, in the 
second and in the third lines of §3§ respectively. We may show that the initial 
monomial of the LHS is no less than the initial monomials of the three terms of 
the RHS with respect to the term order <^_ w w y The initial monomial of the 
LHS biCib k c k dj coincides with that of the 1st term of the RHS. Moreover, all 
monomials appearing in the 2nd and the 3rd term of the RHS have the same 
total degree 5 and (— w, w, 0)-degree 0, and they have a degree at most one with 
respect to bi, bj and b k . This implies that they are less than biCib k c k dj. Thus, 
we have shown that the expression is a standard representation and the S'-pair 
is reduced to zero. 

Since all S-pairs are reduced to zero, the conclusion follows from Buch- 
berger's criterion. Q.E.D. 

Proof of Theorem^ Let I' (a, b, c, d) be the left ideal generated by G' in the ring 
D[a 7 b, c, d]. It follows from Proposition[2]and Theorem 1.2.4] that G' h \ h=1 = 
G' is a Grobner basis of /'(a, 6, c, d) with respect to <(- UJj ^ i o)- Therefore, we 
have 

i"<(-....o) ( 7 '( a > 6 > c ' d )) = < in <<-^,o, (G')> = (&) in D[a, 6, c, d). 

We may replace D[a,b,c,d] by D(a,b,c,d) — C(a,b,c,d)(xij,y k ,r,dij,d k ,d r ). 
Here, C(a, 6, c, d) is the rational function field with variables a = (a pq ),b — 
(h), c = (ci) and d. Then, the following holds: 

in(_ tu , ro )(I'(a, b, c,d)) = (G') in D(a,b,c,d). 

This implies the conclusion. Q.E.D. 
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Proof of the Main theorem^ It follows from ^ that rank (7) > rank (vn/- w W }(I')) 
by Theorem 2.2.1 in [6]. Therefore, we have rank (I) > 2n + 2 by Theorem[2j 
The opposite inequality follows from Theorem 3 of 4 . Q.E.D. 



4 A Proof in the case of n = 1 

In order to clarify ideas of the proof in Section [3J we present a proof of our 
theorem in the case of n = 1. 

The 1-dimensional Fisher-Bingham system of differential equations / C 
C(xi 1 ,xi 2 ,X22,yi,y2,r,dii,di2,d22,d 1 ,d 2 ,d r ) is 

/ = (A n = d n - dl A l2 = d 12 - 9 X 9 2 , 
A22 = d 22 - dl, 
B = dt + dt-r\ 

C12 = xi 2 d\ + 2(x 2 2 - xn)did 2 - xi2&l + y 2 di - y\d 2 , 

E = rd r ~ 2{x n df + x 12 did2 + z 22 <9f ) - (yi<9i + ?/ 2 <9 2 ) - 1). 

The upper bound of rank(J) is 2 • 1 + 2 = 4 as given in 0] Theorem 3]. We 
will show the lower bound of the rank(J) is 4 by using the following general 
inequality 6 S Theorem 2.2.1]: 

rank(J) > in(_ w , w )(/). 

Firstly, we make some change of variables, because it seems to be difficult 
to calculate ia-(- WiW )(I) directly. Let an, ai 2 , a 22 , &i, 6 2 , ci, c 2 , d be constants, 

which will be used as slack variables. We put g — r d ' 'x"" x°^ : 2; 22 2 f where the 
function / is a solution of the system of differential equations /• / = 0. Moreover, 
we make a change of variables 2/1, y 2 by y\ + b\C\, y 2 + 6 2 c 2 respectively. Then, 
the system of differential equations for g is given by I' ■ g = 0, where 

V = (A' n = xndn - x n d 2 - a 3 n , A' 12 = x 12 d 12 - x 12 d 1 d2 - a\ 2 , 
A'22 = x 2 2d 2 2 - x 22 dl - a% 2 , 
B = d 2 + d%-r 2 , 

C[ 2 = x 12 dl + 2(x 2 2 ~ x 11 )d 1 d 2 - x 12 d 2 + {y 2 + b 2 c 2 )di - (y% + b 1 c 1 )d 2 , 
E' = rd r - 2(x 11 d 2 1 + x 12 did 2 + x 22 d%) 

- {{yi + hc^di + (y 2 + b 2 c2)d 2 ) -I- a?). 

We note that rank(J) = rank(I') holds for any set of values of the constants. 

Secondly, we calculate inr WtW ) (I') for the weight vector w = (0, 1, 0, 0, 0, 0) 
where each weight stands for the variables dn, di 2 , 9 22 , d%, d 2 , d r respectively. 
In order to perform Buchberger's algorithm with respect to the (— w, w)-weight 
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order, we need to consider the homogenized system I' h for I': 

I' h = (A* = hxudu - x lx dl - ofi, Ay 2 = hx 12 d 12 - x 12 d 1 d 2 - a?2, 
A 22 = hx 22 d 22 - x 22 d 2 - 022, 
B = df + d\ - rf , 

C 12 = x 12 d 2 + 2(x 22 - x 11 )d 1 d 2 - x 12 d 2 + {y 2 h + b 2 c 2 )dx - (yih + b 1 c 1 )d 2 , 
E' h = hrd r - 2(x 11 df + x 12 d t d 2 + x 22 8 2 ) 

- ((hyi + ha^! + (hy 2 + b 2 c 2 )d 2 ) - h 3 - cP). 

We denote by <^_ w w -, an order in the homogenized Weyl algebra which com- 
pares the total degree first, (— w, w, 0)-degree second, otherwise we apply the 
following block order as a tie breaker: d ^> r ^> {011,012,022} 3> {61 > 
62} > {ci,c 2 } » {2/1,2/2} > d r > {d 12 } » {dn,d 22 } > {di,d 2 } > {2:12} > 
{xii,a;22} 3> h. Here, the symbol > represents the lexicographic order. The 
underlined parts in I are initial terms with respect to <f"_ w w )- They are 
pairwise coprime and their commutators are equal to zero except [B, E' h ] — 
-2hB. From LemmaH we conclude that the set {A% A% A' 2 \, B, E" 1 } is a 
Grobner basis of I' h in D h [an, 012, 022> bi,b 2 , Ci, c 2 , d] with respect to < t ?_ w w y 
In other words, the transformation of dependent and independent variables gives 
us the Grobner basis without adding new elements. Dehomogenizing I' h , we ob- 
tain 

in(_ w ,iu)(/') = (A' n = xudn - x n dl - a\ x , A[ 2 = xi 2 <9i 2 - af 2 , 
A' 22 = x 22 d 22 - x 22 d 2 ~ a\ 2l 
B = d 2 +d 2 -r 2 , 

C' 12 = 2(x 22 - x n )d 1 d 2 + (y 2 + b 2 c 2 )dx - {y x + &ici)<9 2 , 
E' = rd r -2{ Xll di+x 22 dl) 

- {{yi + &iei)<9i + (y 2 + b 2 c 2 )d 2 ) - 1 - d 3 ). 

In this calculation, we regard an,ai 2 ,a 22 ,bi,b 2 ,Ci,c 2 ,d as ring variables with 
the weight 0. As in the proof of Theorem El the equation above holds when 
on, . . . , d are specialized to generic complex numbers. The holonomic rank of 
the (— w, u;)-initial ideal /' := va.(— w>w ) (/') coincides with that of the diago- 
nal system transformed by the same change of variables for I' from /. The 
holonomic rank of I' agrees with that of 

1 = (in = dn -df, Ai2 = 5i2, 
A 22 = d 22 — d 2 , 
B = d[+d 2 2 -r 2 , 

C'12 = 2{x 22 - xu)d 1 d 2 + y 2 di - yid 2 , 

E = rdr- 2( Xll d 2 + x 22 8 2 2 ) - {y x d x + y 2 d 2 ) - 1). 
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Hence, we obtain the following inequality: 



rank(7) = rank(J') > rank(in/_ W)!1 ,)(7')) = rank(J') = rank(I). 
Finally, we show that rank(j) = 4. Proposition [T] tells us that the set 

{All, ^-12) A 2 2, B, C12, E, 

n o/ \fa3 2f} yzdf ~ Vidid 2 - d 2 

D 2 = 2(x 22 - xxx) \d 2 -rd 2 - 



2(x 22 - xn) 



is a Grobner basis of RI with respect to the block order {d r } » {dn > d 22 } 3> 
{d\ > d 2 }, where the tie breaker > represents the graded lexicographic order. 
The set of standard monomials is {1, d\, 82, £?|}. It means that the holonomic 
rank of I is 4. 
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